Abstract. We present results of our diffraction- limited 
mid-infrared imaging of the massive star-forming region 
W3(OH) with SpectroCam-10 on the 5-m Hale telescope 
at wavelengths of 8.8, 11.7, and 17.9 /im. The thermal 
emission from heated dust grains associated with the ul- 
tracompact Hn region W3(OH) is resolved and has a spa- 
tial extent of in the N band. We did not detect the 
hot core source W3(H 2 0) which implies the presence of at 
least 12 mag of extinction at 11.7/zm towards this source. 
These results together with other data were used to con- 
£Sj ■ strain the properties of W3(OH) and W3(H20) and their 
' envelopes by modelling the thermal dust emission. 
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1. Introduction 

The "hot cores" revealed in recent years by molecular 
line investigations are small (< 0.1 pc), very dense (n > 
10 7 cm~ 3 ), and hot (> 100 K) en tities of giant molecular 
clouds (e.g., Cesaroni et al. 1998 ). They are considered to 
be the likely birthplaces of ma ssive stars (e.g., Garay & 



Lizano |1999| , Kurtz et al. 2000 ). Hot cores are frequently 
associated with ultracompact Hn regions (UCHns) which 
are more evolved and accessible to near-infrared (NIR) 
studies of thei r stel lar population (e.g., Feldt et al. 1998 , 
Henning et al. [2001 ). Although the temperatures and sizes 
of hot cores suggest that they might be conspicuous ob- 
jects in the infrared sky, extinction in the mid-infrared 
(MIR, N (10 ^m) and Q (20 ^m) bands) caused by the 
large column densities (A^(H 2 ) > 10 23 cm~ 2 ) must not be 
neglected. 

Conclusions concerning the heating and the stellar con- 
tent of hot cores have to be based on the knowledge of their 
luminosity. This quantity is difficult to estimate since the 
immediate neighbourhood of UCHns often leads to source 
confusion, especially in the far-infrared (FIR) range where 
these objects emit most of their energy and the angular 
resultion of spaceborn observations is as yet poor. Radio 
interferometry at mm/submm wavelengths allow to sepa- 
rate the dust continuum emission of the hot core from the 
free-free radiation of the adjacent UCHn. Such measure- 
ments were used to constrain models of hot cores (Osorio 
et al. 1999 ). High- resolution ground-based MIR observa- 
tions, on the other hand, provide information on the spec- 
tral energy distribution (SED) shortward of the peak flux. 
A corresponding study of the IRc2 source in the Orion 



BN/KL complex (Gezari et al. 1998) illustrates their im- 
portance. More recently, Dc Buizcr et al. (2002) were able 
to detect MIR emission from the hot core G29.96— 0.02 
with a morphology similar to that of the warm ammo- 
nia. We performed high resolution MIR imaging of hot 
cores, including W3(H 2 0) and the neighbouring UCHii 
W3(OH), in order to measure their flux densities or to pro- 



vide at least upper limits. While results for G10.47+0.03 
will be the subject of a forthcoming paper (Pascucci et al., 
in prep.), we present here our findings for W3(H 2 0) and 
the UCHii W3(OH). 

The UCHii W3(OH) is very well-studied in the radio 
domain by continuum and molecular line obser vation s. It 
is located at the distance of 2.2 kpc (Humphreys 1978| ) and 
harbours numerous OH masers. The hot core W3(H2 0), 
also k nown as Turner- Welch object (TW, Turner & Welch 
19841), is situated - 6" 



is situated ~ 6" east of W3(OH). This enigmatic 
source shows an outflow traced by the proper motion of 
H2O masers (Alcolea et al. 1992[ ) and is associated with 
a double-sided radio continuum jet, presumably of syn- 
chrotron nature (Reid et al. 199!:). Recent interferomet- 
ric imaging at 220 GHz by Wyrowski et al. (1999) re- 
vealed another source in the immediate neighbourhood of 
W3(H 2 0), suggesting that the region harbours a cluster of 
protostars. The thermal infrared emission from W3(OH) 



has been studied by Keto et al. ( J1992|) usin g one of the 

Keto et al. 



1987) 
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first MIR array cameras (Arens et al. 
claimed the detection of W3(H 2 0) at the wavelength of 
12.2 ^m with a flux density of 45 ± 10 mJy. W3(H 2 0) and 
W3(OH) were the subject of a continuum and molecular 
line study of van der Tak et al. (2000). Since their results 
are based on single-dish data which do not resolve the two 
objects, they cannot be directly compared to our model 
presented below. 

2. Observations 

The observations were performed using SpectroCam-10, 
the Cornell -built 8-13 /im spectrograph/camera (Hay- 
ward et al. 1993 ) on the 5-m Hale telescope^]. Spectro- 
Cam's detector is a Rockwell 128x128 Si:As BIB array. In 
camera mode, the pixel scale is 0'.'25 with a circular field 
of view of 16". We applied filters with central wavelengths 
and bandwidths (in parentheses) of 8.8 (1.0), 11.7 (1.0), 

1 Observations at the Palomar Observatory were made as 
part of a continuing collaborative agreement between the Cal- 
ifornia Institute of Technology and Cornell University. 
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and 17.9 (0.5) fim. The imaging observations were per- 
formed on 1998 December 27th using the common chop- 
ping/nodding technique with a chopper throw of 20" in 
north-south direction. The wavelet filtering algorithm of 
Pantin & Starck ( 1996| ) was applied to the images. This 
algorithm is useful for recovering faint extended emission 
otherwise hidden in the noise. The high dynamic range of 
the individual frames allowed the application of a shift- 
and-add procedure which includes resampling, yielding a 
final pixel scale of 0"125. To enlarge our field of view, we 
observed two positions at 11. 7 /mi, the first centered on 
W3(OH) and the second offset the east by 10". Lastly, 
the frames were mosaicked. 

The astrometry is based on the radio position of 
W3(OH) which is assumed to coincide with the peak of 
the infrared emission. 

The flux calibration was performed using a Tau as 
reference which was measured at about the same airmass 
(1.2) as the target. For this purpose, the photometric zero 
points from Cohen et al. (1992) and the photometry from 



Cohen et al. (1995) for a Tau were used. The observa- 



tions were performed during photometric conditions, with 
an internal photometric error of less than 2% as estimated 
from the flux variation of the standard star. The 3a sen- 
sitivities (mJy/beam) for the detection of point sources 
in the final images are as follows: 4 (8.8 /im), 6 (11.7/xm), 
and 93 (17.9 /mi). 

3. Results 

The prime results of our observations are contained in 
Fig.|l| which shows the 11. 7 /mi image together with the 
8.4 GH z radio continuum contours from Wilner at al. 
(|l999j beam size 0"2). The two infrared sources in Fig.[j] 
correspond to W3(OH) and the cometary UCHii situ- 
ated north-east of it. A comparison with a Tau (Fig. |] 
insert) shows that W3(OH) is spatially well resolved in 
our diffraction-limited 11.7 /tm image (beam size 0'.'6). The 
lowest contour line of the infrared emission represents the 
3cr detection limit. Clearly, the flux from W3(H2 0) is be- 
low our sensitivity limit. This result is in contradiction 
to the infrared detection of W3(H 2 0) claimed by Keto et 
al. (1992). We note that their conclusion is doubtful since 
the emission they associate with W3(H2 0) peaks only 4" 
instead of 6" east of W3(OH) (see their Fig.l). Further- 
more, their total flux of 45±10mJy contradicts with the 
fact that the emission well exceeds their lowest contour 
level of 100mJy/n". Thus, we conclude that they did not 
detect W3(H20) as well, and were confused by the north- 
eastern cometary UCHn. 

From Fig. [j] it is obvious that there is a good overall 
correspondence between the MIR and the 8.4 GHz radio 
continuum emission. However, there are certain features 
which are different in both maps. The north-eastern trail 
from W3(OH) is more confined in the radio continuum 
and only marginally indicated by the infrared contours. 
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Fig. 1. 11.7/im image of W3(OH) with contours of the 
8.4 GHz radio continuum emission (black lines) from 
Wilner at al. ( 1999 ). The grey contours delineate the 
11.7/mi emission and are spaced by a factor of 2.5 start- 
ing at 3a (lOmJy/D"). W3(H 2 0) is the elongated ra- 
dio source at the reference position. The lower-left insert 
shows the 11.7 /im image of the reference star a Tau using 
a logarithmic brightness scale. 



The peak of the radio emission from the north-eastern 
UCHn is closer towards W3(OH) than its infrared max- 
imum. This is presumably caused by the extinction of 
a dust lane stretching from W3(H 2 0) to the north-west 
which can be seen in the C 17 O(l-0) map of Wyrowski et 
al. ( |1997D . 

For the northeastern UCHn the following fluxes 
(within a 3" aperture) were derived at the three wave- 
lengths: 0.3, 0.6, and 6.1Jy, respectively. 

4. Discussion 



4.1. The UCHn W3(OH) 

The FWHM of the MIR emission from W3(OH) was de- 
rived from Gaussian fits, taking the size of the diffraction- 
limited beam into account. The sizes and flux densities 
of W3(OH) are listed in Tab. [j] The fluxes are based 
on an aperture of 4" diameter. The extent of the MIR 
emission originating from the thermal emission of heated 
dust grains exceeds that of the 8.4 GHz radio continuum 
(FWHM l'/52), indicating that the warm dust is more ex- 
tended than the ionized gas. The variation of the angular 
size in dependence on wavelength can be approximated 
as FWHM(A) ~ A 6±0 - 2 and results from the decline of 
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Table 1. SpectroCam- 
sities of W3(0H) 



10, IRAS-LRS, and MSX flux den- 



\VeLvelengt li 


Peak flux 


± OX ell nux 


T?WP"A/T 
r VV rllVl 


M 


[Jy]/n" 


[Jy] 


["] 


8.8 


2.6 


6.5 


1.96 


11.7 


7.7 


22.4 


2.08 


17.9 


19.5 


93.0 


2.88 


IRAS-LRS 


Total flux 


MSX 


Total flux 


H 


[Jy] 


H 


[Jy] 


8.8 


24.9 


8.28 


n.i 


11.7 


24.1 


12.13 


36.7 


17.9 


89.1 


14.65 


77.1 






21.41 


332.6 



with the high molecular hydrogen column densities of 
1. . . 3.5xl0 24 cm~ 2 inferred from molecular line and con- 



tinuum investigations (Turner & Welch 1984, Wyrowski 
et al. 1997). It suggests that the molecular outflow of 
W3(H 2 0) is very young, i.e. did not fully penetrate the 
hot core yet, and, in addition, might be in the plane of the 
sky. In fact, Fig. 1 from Wyrowski et al. (1999) shows that 
the jet is confined to the region of the hot core. The mod- 
erate expansio n vel ocity of the H2O masers of 20kms _1 
(Alcolea et al. 1992 ) implies a dynamical timescale of only 
500 yr which is consistent with the upper limit on prope r 
motions of the radio jet of 150 kms^ 1 (Wilner et al. 1999 ). 
These velocities are low compared to those of thermal 
radio jets (Anglada 1996 ) and indicate that the outflow 
is presumably slowed-down by the high-density environ- 
ment. 



the temperature with increasing distance from the heat- 
ing star(s). 

The comparison of our flux densities of W3(OH) with 
other estimates allows conclusions on the influence of dif- 
ferent beam sizes. For this purpose, we retrieved the IRAS- 
LRS spectrum, identified W3(O H) in the MSX point 
source catalog (Egan et al. 1999 ), and retrieved an ISO- 
LWS spectrum from the data archive. 

The LRS spectrum was integrated according to the ap- 
plied passbands. It is obvious that the 8.8 /jm fluxes given 
in Tab.|l| considerably exceed our value. This can be ex- 
plained by ubiquitous emission attributed to Polycyclic 
Aromatic Hydrocarbons (PAHs) surrounding the UCH11 
which strongly contributes to the flux in the large aper- 
tures of IRAS and MSX. Pronounced 7.7 and 8.6 /jm PAH 
bands can be misleading in ground-based derivations of 
the optical depth of the 9.7 silicate feature (e.g., Roelf- 
sema et al. 1996). 



4.2. The hot core W3(H 2 0) 

Our attempt to detect W3(H20) in the IR was stimu- 
lated by the presence of an outflow. Generally, IR emis- 
sion can escape in outflow lobes primarily due to scattering 



e.g., Fischer et al. 1996). An example is NGC6334 I(N), 



a presumed high-mass Class object, for which Sandell 
(Sandell |2001 ) rendered the detection of IR emission im- 
possible because of the extremely high extinction derived 
from mm/submm maps. However, this sourc e is asso- 
ciated with NIR emission (Tapia et al. 1996 , Megeath 
& Tieftrunk |l999[ ) obviously originating from the blue- 
shifted lobe of its outflow. 

The flux densities from W3(H20) in the absence of any 
intervening absorbing matter can be estimated from the 



temperature map given by Wyrowski et al. (1997). The 



expected peak surface brightness amounts to 2170 Jy/d" 
at 11.7/zm. Together with our 3<7 sensitivity, this yields 
a lower limit to the extinction at this wavelength of 
12 mag. Our failure to detect this source is consistent 



5. Modelling the dust emission 

One of the questions concerning W3(OH) and W3(H 2 0) 
is related to their individual luminosities which directly 
translates to the nature of the internal heating sources. Al- 
though it is reasonable to assume that the internal heating 
is due to stars of intermediate or high mass, no direct con- 
firmation of their presence exists. The luminosity might 
also, at least partly, result from accretion. The high col- 
umn densities toward W3(H 2 0) render it very difficult to 
figure out whether the mass accretion is due to infall from 
a circumstellar disk or stellar mergers like in the scenario 
of Bonnell et al. ( 1998] ). The high optical depths will lead 
to a thermalization of the released energy irrespective of 
the accretion mechanism. 

The derivation of the individual luminosities has to 
be based on the decomposition of the SED. While our 
measurements provide constraints for the MIR, the lack 
of spatial resolution in the FIR does not permit to sepa- 
rate both components. An upper limit of ^1800 Jy on the 
50 ^m flux o f W3(H 2 Q) has been established by Campbell 
et al. (1989) from KAO scans. In the mm/submm range, 
aperture synthesis measurements show that the 1.3 mm 
emission from W3(H 2 0) is due to dust radiation while 
free- free emission dominates the flux of W3(OH) at this 
wavelength (Wyrowski et al. |l999j ). The ISO-LWS spec- 
trum between 42 and 140 /im can well be fitted by a mod- 



ified black-body (35 K and k(v) 



'). The association 



of W3(H 2 0) with cold dust suggests that the MIR/FIR 
excess in the ISO spectrum for A > 140 /im is presumably 
due to the hot core. 

We modelled the thermal dust radiation from W3(OH) 
and W3(H 2 0) in order to figure out how much the hot 
core contributes to the FIR excess seen in the ISO-LWS 
spectrum and which sensitivities are required for its de- 
tection in the MIR. The input SED is based on fluxes 
coming from this work, from MSX point source catalogue 
at 21.4 /im (see Tab. 0), from the IRAS -LSR and ISO- 
LWS spectra, and from Wyrowski et al. 1997 as well as 
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Table 2. Radiative transfer model parameters 

W3(OH) 
W3(H 2 0) 



Source 


L* 




Rout 


M d 




[10 4 L Q ] 


[AU] 


[AU] 


[M ] 


W3(OH) 


8 1 


2270 2 ' 3 


57000 4 


5.0 


W3(H 2 0) 


2 4 5,6 


190 


24000 5 


3.0 7 



Ref: (1) Wink et al. 1994, (2) Chini et al. 1986, 
(3) Campbell et al. 1989, (4) Zeng et al. 1984, 
(5) Osorio et al. 1999, (6) Thompson 1984, 
(7) Wyrowski et al. 1997 



Wilner et al. 1995 for the mm wavelenghts. To solve the 
radiative tra nsfer problem we used the ID-code of Manskc 
& Henning ( 1998 ) assuming spherically symmetric shells 
like in previous works. As for the dust composition, we use 
graphite, silicate, and iron with optical consta nts ta ken 
from Dorschner et al. ( 1995 ) and Draine & Lee ( 1984 ) to- 
gether with a MRN-type size distribution (Mathis et al. 
1977 ). The model parameters and the related references 
are given in Tab. ^: here L» is the stellar luminosity, Rj„ 
and R ou t are the inner and the outer radius of the en- 
velope, and Md is the dust mass. We underline that the 
inner radius of the dust shell of W3(OH) coincides with 
the outer radius of the UCHii while for the hot core, Ri„ 
is derived from the dust sublimation temperature (~ 1000 
K). T o model the W3(H 2 0) SED we follow Osorio et al. 
(1999) in adding an accretion luminosity of 4.2 x 10 4 L 
to the stellar luminosity and a free-fall density shell. For 
W3(OH) the density profile that better fits our observa- 
tions is gaussian and yields an almost constant density 
distribution till ~ 20000 AU. The resulting opacities at 
the mm wavelenghts varies as v 15 . 

Fig. ^ shows the individual SEDs as well as their superpo- 
sition. Our model produces the following FWHMs at the 
observed wavelenghts: 1'.'6 at 8.8 /mi, 2'.'0 at 11.7/im and 
3'.'1 at 17.9 /im. They are quite in agreement with those 
given in Tab. [l] which supports the view that the parame- 
ters of the radiative transfer model are representative for 
the actual conditions. From the SED of W3(H 2 0) we con- 
clude that deep observations in the Q band might be able 
to detect the MIR emission from the hot core. 



6. Summary 

1. The thermal infrared counterpart to the hot core 
source W3(H2 0) was not detected at any wavelength 
of our observations. This revises the finding of Keto et 
al. (1992) which most probably resulted from confus- 



ing W3(H 2 0) with the UCHii northeast of W3(OH). 
We derived a lower limit for the extinction at 11.7/im 
towards W3(H 2 0) of 12 mag. 

Our diffraction-limited imaging resolved the thermal 
emission from W3(OH) and clearly indicates a wave- 
length dependence of its apparent size. 
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Fig. 2. SEDs of the individual sources (W3(OH) - dot- 
ted line, W3(H20) - dashed line) and of their sum (solid 
line) . Asterisks with error bars indicate the measured flux 
densities. Diamonds represent the model fluxes for the ob- 
servational beam sizes. Arrows mark the detection limits 
for W3(H 2 0). 

3. The comparison of our 8.8 /tm flux to those measured 
with IRAS and MSX led to the conclusions that PAHs 
surround W3(OH). 

4. In accordance with our thermal infrared imaging 
and mm/submm studies, the ISO-LWS spectrum of 
W3(OH) might be decomposed in two components 
with W3(OH) being the hotter, more evolved, object 
while W3(H 2 0) dominating at mid/far- infrared wave- 
lengths. 
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